Habitat suitability and connectivity modeling predict genetic population structure and priority control areas for invasive nutria (Myocastor coypus) in a temperate river basin

The nutria (Myocastor coypus), also known as the coypu, is a semi-aquatic, invasive rodent native to South America that causes damage to natural riverine and wetland habitats in many parts of the world, including South Korea. Understanding habitat use, connectivity, and gene flow of nutria populations is critical for the sound management of local and regional ecosystems. Here, we assessed habitat suitability and connectivity in relation to the genetic structure of nutria populations in the Nakdong River Basin of South Korea. A total of 321 nutria occurrence sites and seven environmental variables were used to perform ensemble habitat suitability modeling using five species distribution models (SDMs), including boosted regression trees, maximum entropy model, random forest, generalized linear model, and multivariate adaptive regression splines. Using graph and circuit theory approaches, we assessed the population gene flow and current flow betweenness centrality (CFBC) of suitable habitats derived from the ensemble SDM. All SDMs performed well with a range of test AUC values from 0.962 to 0.970 (mean = 0.966) with true skill statistic values over 0.8. The minimum temperature of the coldest month, mean temperature of the warmest quarter, precipitation of the driest quarter, and distance from water bodies were important predictors in nutria habitat modeling. Nutria population gene flow was significantly correlated with the least-cost path distance on a cost resistance surface based on ensemble habitat suitability modeling and roads (Mantel’s r = 0.60, p < 0.05). Finally, the CFBC positively correlated with the genetic diversity of nutria populations was used to identify priority control areas. Habitat suitability and connectivity modeling not only revealed environmental conditions and areas that support the survival and spread of nutrias, but also improved our understanding of the animals’ genetic population structure, thereby indicating priority areas to target for eradication.


Introduction
Invasive alien species (IAS) are one of the largest threats to native biodiversity, along with habitat loss and climate change [1][2][3]. IAS compete with other native species for limited resources, cause disease and extinction of indigenous flora and fauna, and lead to degradation and displacement of native habitats [4][5][6]. This can induce significant ecosystem disruption and socio-economic damage by altering ecosystem services or impacting infrastructure [7,8]. Moreover, the management of IAS, either through control or eradication, usually requires considerable labor and financial resources [9].
Understanding the factors driving invasion success, distribution, spread pattern, and impacts is critical for the effective and efficient management of IAS [10]. Although preventing new IAS from arriving is the most cost-effective solution, the introduction of alien species and wildlife trading are continuously increasing, facilitating the spread of IAS [11,12]. Once an IAS is established and has spread, identifying current and potential habitats and their connectivity becomes an essential task to effectively control their impact in ecosystems [13]. Because the spread of IAS is constantly increasing worldwide, despite global efforts to limit biological invasions [14], prioritizing high-risk IAS based on their potential ecological and socio-economic impacts is crucial.
The nutria (Myocastor coypus), also called the coypu, is a large, semi-aquatic rodent native to temperate and subtropical South America. The species was initially introduced to Africa, Asia, Europe, and North America for fur and meat farming. Due to declining demand for nutria products, nutrias have been released and flourished in natural riverine and wetland habitats in Korea and other parts of the world. The nutria has generalist feeding habits, mostly fond of aquatic plants and roots, thus damaging native vegetation and crops, resulting in changes in soil structure, and turning valuable wetlands into open water [15,16]. Nutria burrowing and foraging behaviors further aggravate wetland and biodiversity loss by accelerating soil erosion and competing with native species, such as the muskrat in North America [17][18][19].
Because the nutria is recognized as a pervasive ecological threat, its establishment and spread has been identified as one of the top priority environmental issues in many places [20][21][22]. Accordingly, the nutria has also been listed among the top 100 of the world's worst IAS [23]. However, the detailed distribution pattern and, by extension, the landscape genetics of the nutria remain unclear, which may provide vital clues for effectively detecting and controlling its spread. Some studies have related the occurrence of nutrias with environmental factors such as climate and topographic conditions to predict habitat suitability [24][25][26] or have investigated the distributional shifts of nutria populations [20,26,27]. Other studies have characterized the structure and genetic diversity of nutria populations for management purposes [28,29]. However, no comprehensive study has yet examined the gene flow and habitat connectivity of nutria populations as well as their habitat distribution, which can serve as a key functional element in the development of spatial management plans and strategies for nutria control. Connectivity and gene flow, in turn, shape the genetic structure pattern of a species, which is of profound evolutionary and ecological importance [30,31]. Overall, the habitat suitability and gene flow network for nutrias are poorly understood, which hinders the implementation of IAS management to effectively isolate and eradicate nutria populations.
To overcome these limitations, we developed nutria habitat suitability and connectivity models based on the occurrence of nutrias, bioclimatic and environmental data, and population genetic data. Focusing on this invasive species, we demonstrated how models combining habitat and genetic data can be used to predict species distribution and spatial population structure. We specifically aimed to assess habitat suitability, gene flow connectivity, and the habitat network of nutria populations in the Nakdong River Basin of South Korea. Using ensemble species distribution models along with graph and circuit theory approaches, we identified important environmental factors influencing nutria distribution, assessed habitat suitability and gene flow connectivity, examined the functional habitat network, and prioritized key areas to control the spread of this IAS in South Korea. The overall research process and outcomes of this study are shown in Fig 1.

Study area
The study area was comprised of the Nakdong River Basin and surrounding islands in Korea, lying between 34.5-37.3˚N, and 127.4-129.6˚E, WGS 84 (Fig 2A). The Nakdong River is the longest river in South Korea (approximately 510 km) ( Fig 2B). The study area was approximately 31,669 km 2 with an average elevation of approximately 266 m above mean sea level. Forests cover 68.4% of the Nakdong River Basin; upland fields and paddy fields cover 7.2% and 10.4% of the basin, respectively; with the remaining areas covered by grasslands, bare lands, and urban areas (0.1-6.9%) [32]. The study area had a typical temperate climate with four distinct seasons affected by East Asian monsoons.

Nutria occurrence sites
Between 2016 and 2020, field surveys of the study area were conducted at a walking speed of 0.5-2.0 km/h, more than twice every quarter, i.e., each season, by 10-15 wildlife experts of the Nakdong River Basin Environmental Office (Nakdong River BEO, Ministry of Environment of Korea). The locations of nutria occurrences based on direct and indirect observations, such as dens, droppings, and footprints, were noted during the field surveys. All sites were recorded using GPS coordinates. To minimize location overlaps and maintain compatibility with spatial data resolution, all nutria occurrence sites were separated by a minimum distance of 1 km. A total of 321 nutria occurrence sites were obtained ( Fig 2B).

Species distribution modeling
Seven bioclimatic and environmental variables were selected based on their relevance to nutria habitat distribution according to previous studies [25][26][27]33] (Table 1). Among them, three bioclimatic variables with a spatial resolution of 30 arc-seconds (1 km) were included: the . The genetic data location coordinates were obtained from the study by Kim et al. [28]. River and stream data were derived from the OpenStreetMap database (http://www.openstreetmap.org).
https://doi.org/10.1371/journal.pone.0279082.g002 minimum temperature of the coldest month, the mean temperature of the warmest quarter, and the precipitation of the driest quarter, averaged from 1970-2000 [34]. As a ground vegetation biomass variable, the median normalized difference vegetation index (NDVI) at 30 m spatial resolution was generated using all Landsat-8 operational land imager (OLI) data from 2016-2020 via the Google Earth Engine (GEE) platform [35]. Distance from water bodies was included as a predictor of foraging activity [36]. Two anthropogenic factors were also used as human impact proxies, including distance to roads and human population density. The Open-StreetMap (OSM) database (as of November 2017) provided road vector data, including motorways, trunks, and primary, secondary, and tertiary roads [37]. As a population density variable, population grid data of 1 km resolution was obtained from the geospatial information platform (http://map.ngii.go.kr) for October 2018. To maintain consistency with other data, the NDVI variable was resampled to a 1 km resolution using bilinear interpolation. To examine multicollinearity, Pearson's cross-correlations were calculated between all variable pairs based on the 321 nutria occurrence sites and 10,000 randomly generated points across the study area. No variables were highly correlated (|r| < 0.7 for all pairs) [38]. An ensemble of predictions of five species distribution model (SDM) algorithms was employed to identify the major environmental factors influencing nutria distribution and potential suitable habitats, using Software for Assisted Habitat Modeling (SAHM) v.2.0.1 for VisTrails v.2.2.3 [39]. The SAHM software provides a workflow interface that enables easily to fit SDMs, to calculate model performance metrics, to explore variable importance and response curves, and to perform model selection. The five SDM models included boosted regression trees (BRT), maximum entropy model (MaxEnt), random forest (RF), generalized linear model (GLM), and multivariate adaptive regression splines (MARS). A presence-only method using background data was employed to fit the five models [40].
The 321 nutria occurrence sites were combined with environmental factors in the SDMs for habitat suitability predictions. A total of 10,000 background data points were randomly generated, excluding nutria occurrence sites within the study area. For BRT, alpha = 1, bag fraction = 0.75, number of folds = 3, and default settings for the other values were used. For MaxEnt, the default settings in the SAHM package were used. For RF, the out-of-bag prediction error was minimized using the tuneRF function. For GLM, a bidirectional stepwise procedure based on the Akaike information criterion (AIC) was used, considering all interactions and squared terms. For MARS, Mars degree (Friedman's μ) = 1 and generalized cross-validation penalty per knot = 2.0 were used [41].
All SDMs were evaluated based on a 70-30% train-test data split using the area under the receiver operating characteristic (ROC) curve (AUC) and true skill statistic (TSS) [42]. The AUC (a threshold-independent evaluation metric) ranges from 0 to 1, with higher values indicating better model performance [43]. The TSS (a threshold-dependent evaluation metric) ranges from -1 to 1, with higher values indicating better model performance [42]. The TSS was calculated using the threshold that maximized the sum of sensitivity and specificity in the ROC curve [44,45]. Ten-fold cross-validation was applied to the training dataset for model selection. Only models with an average test AUC value > 0.75 and TSS value > 0.4 were selected for further analyses [46].
The variable importance for each model was examined using the change in AUC values with and without the independent variable, but using all other variables [47]. The increase in AUC (ΔAUC) was calculated when each variable was permuted for each model. Variable importance was then ranked by the mean ΔAUC of the five models and model response curves were generated for the four most important variables with ΔAUC � 0.05.
To obtain a robust ensemble model estimate, model predictions were consolidated using the committee averaging method [48] in which model predictions for occurrence probability were translated into binary form (1 = suitable, 0 = unsuitable) using the threshold that maximized the sum of sensitivity and specificity [44,45]. The binary prediction models were then averaged to create a final ensemble model of habitat suitability [49].

Gene flow and isolation modeling
The fixation index (F ST ) is a widely used genetic distance measure for genetic differentiation [50] which occurs when gene flow between populations is restricted. Pairwise F ST values of eight nutria populations in the Nakdong River Basin of South Korea, comprised of 93 individuals each typed at 12 microsatellite loci, were obtained from the study by Kim et al. [28] ( Fig  2B). These authors reported that variation in F ST values was not explained by pairwise Euclidean distance (i.e., isolation-by-distance, IBD). Therefore, to estimate the presence of genetic isolation by other advanced distance measures, the F ST genetic distance was evaluated against the least-cost distance and the resistance distance among population pairs [51,52]. The leastcost distance (LCD) is the accumulated cost distance of the optimal pathway connecting population pairs. The resistance distance (i.e., isolation-by-resistance, IBR) considers all possible pathways connecting population pairs based on circuit theory [52]. To measure the two distance metrics, a 50 m resolution resistance surface was generated by assigning exponentially decreasing cost values (ranging from 1-100) to the habitat suitability values from the ensemble model as follows: a cost value of 100 assigned to 0 (least suitable), 40 to 0.2, 15 to 0.4, 5 to 0.6, 2 to 0.8, and 1 to 1 (most suitable). In addition, a cost of 10,000 was assigned for major roads since these seem to be one of the most disruptive barriers to nutria migration [53]. The effect of roads on gene flow was also examined by either including or not including the road cost.
LCDs were calculated using PATHMATRIX software [54]. Resistance distances were measured using Circuitscape v4.0.5 software [55]. The significance of correlations between F ST and distance measures were evaluated using Mantel tests [56], and were computed using the vegan package in R with 10,000 permutations [57].

Habitat network analysis and area prioritization
A habitat network map was created to identify priority areas for nutria control. Graph and circuit theory approaches were used to map the habitat network, providing effective methods to evaluate multiple aspects of habitat connectivity and priority areas [58,59]. The habitat network (i.e., graph) was defined as a set of suitable habitat patches (nodes) and connections between them (links). In the final ensemble SDM, nutria occurrence probabilities above 0.5 were considered suitable (at least three of the five SDMs predicted as suitable). We rasterized the suitable pixel center points of the ensemble prediction map at a 50 m resolution, and defined them as habitat patches.
As a precautionary principle, the degree of connectivity should be overestimated rather than underestimated when connectivity is viewed negatively [60]. Thus, least-cost modeling to find the optimal connectivity solution is considered a better way to measure connectivity for invasive species. For this reason, links in the network, representing the functional relationships between habitat patches in terms of potential migration or movement of nutrias, were defined by a cumulative cost distance using the 50 m resolution cost surface, which used also for the gene flow and isolation modeling. The cumulative least-cost distances between habitats were within the maximum cost distance of 250,000 for computational efficiency. The habitat network was assumed to be a thresholded complete graph that directly connected all pairs of nodes with less than a specified threshold cost distance. Nutria dispersal and home range distances, also known as the population's spreading rate, are known in Euclidean form, ranging from hundreds of meters to tens of kilometers [61]. Moreover, the home range size and dispersal distance of nutrias can be largely dependent on resource availability and environmental conditions [61,62]. Thus, multiple cost distances from 5,000-25,000 with intervals of 5,000 were employed to threshold the habitat network.
The connectivity metric (current flow betweenness centrality (CFBC) [63]) was calculated to determine priority areas for nutria control. CFBC applies a circuit theory algorithm [64]. Links on a regular lattice (e.g., raster grid) correspond to resistors, each of which has its own resistance. In this study, current (i.e., flow of charge), which can be measured through a resistor, represented the movement probability of a random walker through least-cost paths. CFBC is determined only by the resistance and number of links. The metric requires only one parameter (e.g., patch capacity) that influences the distribution of CFBC along the nodes and links of the graph. In this study, the capacity acted as the amperage of the source patches. According to the patch definition, the amperage was set to 1 for all patches (β = 0). The CFBC estimates of each habitat network thresholded at the different cost distances were measured at the patch level. To determine the optimal cost distance scale in relation to the genetic structure of nutria populations, Pearson correlations were computed between expected heterozygosity (He, also called gene diversity) and CFBC measured at each threshold cost distance. Due to lack of normality, the CFBC estimates were log-transformed for Pearson correlation analysis. He values of eight nutria populations were obtained from the study by Kim et al. [28]. Each nutria population's CFBC corresponded to that of the nearest habitat patch to each genetic data locality. High-priority control areas were identified as those within the top 5% of CFBC derived at the most significant scale relating nutria genetic diversity with habitat connectivity. Habitat network analyses were conducted using Graphab v2.6.1 software [65]. ArcGIS v10.4.1 software was used for mapping [66] and IBM SPSS Statistics 25 [67] to perform correlation tests.

Results
All five SDMs performed well in predicting nutria occurrence, with test AUC values ranging from 0.962-0.970 with a mean of 0.966. MaxEnt (AUC = 0.970) had the highest predictive capability, followed by RF (AUC = 0.968), BRT (AUC = 0.965), MARS (AUC = 0.963), and GLM (AUC = 0.962). The TSS values of all models ranged from 0.825-0.888. When assessing model performance using a ten-fold cross-validation on the training dataset, all five SDMs had good predictive power, with a mean test AUC � 0.954 and TSS � 0.442. Thus, the predictions of all five SDMs were combined into an ensemble SDM using the committee averaging method.
Among the seven predictors, the minimum temperature of the coldest month was the most important variable (Table 2). For most SDMs, this predictor generally had a positive influence on nutria occurrence at minimum temperatures higher than approximately -6˚C, while the probability of occurrence was nearly 0 at minimum temperatures lower than approximately -6˚C (Fig 3). As the second most important variable (Table 2), the mean temperature of the warmest quarter generally had a positive effect on nutria occurrence at mean temperatures higher than approximately 23.5˚C (Fig 4). The precipitation of the driest quarter and distance from water bodies were identified as other important variables ( Table 2) negatively correlated with nutria occurrence (Figs 5 and 6). Median NDVI, distance to roads, and population density were not important predictors of nutria occurrence (i.e., ΔAUC values < 0.01). Fig 7 presents the ensemble map derived from the five SDMs for nutria occurrence in the Nakdong River Basin, Korea. Suitable nutria habitats predicted by all and � 3 models covered 2.84 and 9.73% of the total basin area, respectively. Conversely, unsuitable habitats predicted by all and � 3 models covered 83.85 and 90.27% of the total area, respectively. The most suitable nutria habitats (i.e., core habitats) were in the middle and downstream lowland areas of the Nakdong River Basin (Fig 7).
The F ST genetic distance was significantly correlated with LCD based on habitat suitability and road cost values (Mantel r = 0.600, p < 0.05), but not with the IBR (Mantel r = 0.154, p = 0.333) (Fig 8). Least cost pathways connecting nutria populations mainly followed the Nakdong River and its major tributaries (Fig 9). However, genetic distance was not correlated with LCD when road cost was not included (Mantel r = 0.110, p = 0.359). The He values of the eight nutria populations were significantly correlated (p < 0.05) with their habitat CFBC measured at distance thresholds of 5,000-15,000 cost units (Fig 10A). The most significant correlation between the He values of the eight nutria populations and CFBC (r = 0.758, p < 0.05) occurred when the distance threshold was set at 5,000 cost units (Fig 10B). Thus, the distance threshold for the network analyses was set at 5,000 cost units. The map of the CFBC pinpointed high-priority control areas, including the top 10% of habitat nodes with high connectivity along the midstream areas (Fig 11). Riparian habitats in the mainstream section of the Nakdong River and its major adjacent tributaries, approximately 100 km long, were identified as key control areas for the prevention of nutria spread and eradication.

Discussion
This study is the first comprehensive modeling assessment to predict habitat distribution and priority control areas for invasive nutrias in a temperate river basin. The study findings indicate that habitat suitability was largely determined by the minimum temperature of the coldest month, mean temperature of the warmest quarter, precipitation of the driest quarter, and distance from water bodies. Further, the level of least-cost connectivity based on an ensemble model and road network was significantly related to the genetic differentiation of nutria populations. In addition, the genetic diversity of nutria populations was positively related to the degree of node connectivity in the species habitat network. Importance analysis of environmental variables highlighted two temperature-related climatic factors-the minimum temperature of the coldest month and the mean temperature of the warmest quarter-as being crucial in shaping habitat suitability for the nutria. Nutria occurrence was more likely in areas of mild climate with the minimum temperature of the coldest month and the mean temperature of the warmest quarter above approximately -6 and 23.5˚C, respectively. This finding concurs with that of other studies reporting temperature as a potential limiting factor of spread [20,26] since nutrias are particularly vulnerable to cold winters [27]. Indeed, nutrias tend to thrive and propagate rapidly in areas with mild winters [62,68]. Although cold weather has been identified as a significant factor limiting the reproduction and range expansion of nutrias [25,69,70], ongoing climate change may lead to range expansion of the species [26]. The precipitation of the driest quarter and distance from water bodies also played an important role in shaping the distribution of the species. The probability of nutria occurrence decreased with increasing levels of precipitation in the driest quarter. In the study area, the driest quarter broadly coincides with the cold season. Thus, the precipitation of the driest quarter can be considered a proxy for snow cover or depth, which may hinder animal movement and feeding [71]. Distance from water bodies also had a negative effect on nutria distribution. This result is consistent with the ecological characteristics of nutrias, which rely on riparian habitats. Previous studies have also indicated the importance of water resources (e.g., rivers and streams) for the distribution of this species [25,27,33]. Meanwhile, the NDVI variable was not important in the SDMs, even though vegetation cover provides food and shelter [33,61,72]. Considering the ecological characteristics of the nutria, the combined effect of vegetation cover and water resources on habitat suitability should be evaluated using methods such as the ecological niche factor analysis (ENFA) [33]. The other variables, including distance to roads and population density, also had lower importance values, as reported in other studies [26,33].
The ensemble distribution map (Fig 7) revealed that nearly 10% of the Nakdong River Basin in Korea can be considered suitable nutria habitats, predicted by over half of the SDMs. The map also illustrated that nutrias prefer lowland riverine areas compared to mountainous areas [33]. However, the predicted suitable habitats were confined to midstream and downstream riverine areas of the basin with relatively mild winters. Due to continuous global warming, temperatures in the basin are expected to rise, most notably in upstream regions in winter [73]. Although no occurrence was recorded in the upstream regions in this study or in other studies [20,27], upstream areas adjacent to suitable nutria habitats are at risk of invasion. Therefore, these areas warrant further attention and monitoring.
Population structure in landscape genetics is commonly assessed by IBD, examining the relationship between Euclidean inter-population geographic distance and genetic distance   https://doi.org/10.1371/journal.pone.0279082.g006 [74]. However, the current study applied LCD and IBR models as more sophisticated landscape genetic approaches to evaluate the effects of spatial habitat suitability and roads on genetic connectivity. Although Kim et al. [28] reported that variation in F ST values for eight nutria populations was not explained by IBD, a significant positive correlation was determined in the current study between genetic and LCD on a resistance surface based on an ensemble habitat suitability model and roads (Fig 8A). This result indicates that gene flow of nutria populations is highly constrained by the spatial configuration of habitat suitability and roads. Further, least cost path modeling suggested that potential gene flow routes mainly follow the Nakdong River and its major tributaries. In addition, genetic distance was not significantly correlated with LCD when major roads were not considered movement barriers, whereas distance to roads had a low importance in SDMs. Overall, these suggest that nutria can occur close to roads, but roads are significant barriers to nutria migration. Thus, suitable habitat areas with less road coverage should be prioritized for the control and prevention of nutria spread.
Circuit theory-based IBR often predicts population genetic structures better than IBD and LCD [75]. Circuit theory, which is based on random-walk theory, assumes that animals have no knowledge of the landscape (i.e., the relative resistance of landscape features) beyond their immediate surroundings [76]. In contrast, least-cost modeling assumes that animals have knowledge of the landscape and associated costs in order to choose the lowest cost path [51]. Interestingly, unlike LCD, IBR was not supported in this study (Mantel r = 0.154, p = 0.333). Nutrias spread or disperse mainly along rivers and tributaries [27] and tend to remain in the vicinity of their natal area for their entire lives (i.e., philopatric behavior) [77]. These characteristics suggest that a cost resistance surface mainly based on habitat suitability ranging from 0 to 1 may not be appropriate for IBR to evaluate the contribution of multiple pathways to the gene flow of this species.
High habitat connectivity may increase genetic diversity through increased gene flow [78]. Accordingly, the study results demonstrated that the degree of habitat connectivity based on network flow (i.e., CFBC) was positively correlated with the genetic diversity of nutria populations. Moreover, habitat connectivity was related to the genetic diversity of the species in a scale-dependent manner. He values of eight nutria populations were significantly correlated (p < 0.05) with their habitat CFBC, measured at distance thresholds of 5,000 to 15,000 cost units (Fig 10A). This may correspond to the range of scales at which nutrias perceive habitats that are functionally connected in the middle and downstream regions of the Nakdong River Basin, which might be associated with species traits, such as movement behavior and perceptual range, as well as environmental conditions. Despite their philopatric behavior, nutrias are known to travel tens of kilometers when resources are limited [25,79]. However, due to suitable environmental conditions for nutria survival in this region, the most significant correlation was found at a relatively short cost distance of 5,000 units (corresponding to a Euclidean distance of approximately 2.5 km).
Priority areas to control the establishment and spread of nutrias were explicitly identified using a CFBC map as the backbone of habitat connectivity in the Nakdong River Basin (Fig  11). Key control areas of high habitat connectivity included riparian habitats in the mainstream section of the Nakdong River and its major adjacent tributaries. These areas pinpointed where surveillance and control efforts, such as hunting and trapping, should be targeted to effectively and efficiently eradicate and isolate nutria populations. Surveillance and control of these areas, especially those close to high biodiversity and protected areas, are critically important for ecosystem conservation and restoration.

Conclusions
The nutria poses a serious ecosystem threat due to its impacts on biodiversity and water resources. Targeted actions against species invasion should be based on a reliable, integrated, and spatially explicit understanding of habitat requirements and ecological processes, as well as available control options. To this end, this study combined ensemble habitat suitability modeling, landscape genetics, and habitat network analysis to create a functional habitat network map as an integrated approach for the prevention of nutria spread and its eradication. This habitat network map is expected to play a pivotal role in identifying priority control areas Least-cost paths connecting pairs of eight nutria populations in the Nakdong River Basin (base map data source: https://gadm.org/). The middle and downstream area was enlarged. The genetic data location coordinates were obtained from the study by Kim et al. [28]. River and stream data were derived from the OpenStreetMap database (http://www.openstreetmap.org). https://doi.org/10.1371/journal.pone.0279082.g009

PLOS ONE
and guiding the planning and management of invaded and adjacent non-invaded riparian areas to reduce the distribution and spread of nutrias. This study provides a cost-effective and efficient approach to better allocate resources for continuous surveillance, monitoring, and control of invasive species at local or regional levels. Because species distribution variability  habitat nodes measured at a cost distance threshold of 5,000 within the middle and downstream regions of the Nakdong River Basin. Circles represent habitat nodes. Red and yellow circles indicate priority habitat nodes of high connectivity for invasive species management. (B) Enlarged view of a key control area including high-priority habitat nodes along the Nakdong River and its major adjacent tributaries. Links connecting habitat nodes are also shown. River and stream data were derived from the OpenStreetMap database (http://www.openstreetmap.org).
https://doi.org/10.1371/journal.pone.0279082.g011 may occur at relatively smaller scales, in the future, high resolution local environmental data are required to perform more accurate habitat suitability and connectivity modeling for nutria control. Further research should also be directed to better understand individual movement and space-use behavior patterns of this species across spatiotemporal scales (i.e., movement ecology) in a changing landscape and climate.
Supporting information S1 File. The WGS84 geographic coordinates of nutria occurrence locations used to fit the five species distribution models. (ZIP)